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THEORY OF THE ANISOTROPIC HEISENBERG FERROMAGNET 

by Lawrence Flax 
Lewis Research Center 

SUMMARY 

The development of a formalism for the investigation of phase transitions of the 
anisotropic Heisenberg ferromagnet is presented. By the use of Green's function theory 
the magnetization is determined by several decoupling schemes. A method of calculating 
the generalized Watson sums is also formulated. 

The thermodynamic quantity to be studied is the magnetization at and near the Curie 
point. By use of various decoupling schemes, different phase transitions are generated. 
An analytical formulation is developed to determine at what temperature the transition 
takes place. 


INTRODUCTION 

As the temperature of a ferromagnet is increased, the magnetization will decrease 
until, at a temperature known as the Curie temperature T £ the sample becomes para- 
magnetic. The process by which this change from ferromagnet to paramagnet occurs is 
referred to as the ferromagnetic phase transition. The fundamental problem in the 
theory of ferromagnetism is the theoretical description of the phase transition. The 
theory of this phase transition, as well as many other types of phase transitions occur- 
ring in nature, has been the subject of many investigations. These investigations have 
shown that electrical, mechanical, and many thermodynamic properties of a material 
are altered when the material undergoes a phase transition. 

One way of describing the ferromagnetic phase transition is to look upon it as an 
order-disorder transition in the material. The property of order in a system describes 
the correlation between the properties of two atoms, even when the atoms are widely 
separated in space. The spontaneous magnetization present in a ferromagnet is due to 
the correlations between the spins of the atoms. Above the Curie temperature the long- 
range order is gone, and the system is paramagnetic. The description of this change 
from an ordered to a disordered state is the description of the thermodynamic phase 
transition. 



Thermodynamic phase transitions can be classified by the discontinuities that occur 
in the Gibbs thermodynamic potential. From an experimental standpoint, it is not easy 
to decide what types of transition are present because of the rapid variation of a thermo- 
dynamic variable such as specific heat or entropy, near the critical temperature. It is 
often difficult to distinguish, for example, whether a transition in specific heat against 
temperature has a finite, infinite, or cusp discontinuity. 

This report is concerned with the investigation of the thermodynamic properties of 
the anisotropic Heisenberg ferromagnet. The thermodynamic quantity studied was the 
magnetization at and near the Curie point. The statistical method used to determine the 
thermodynamics of the system was the double time thermodynamic Green's functions. 
The objectives of this report are threefold: One is to investigate the various decoupling 
schemes for the Green's function theory of the anisotropic Heisenberg model. The sec- 
ond purpose is the evaluation of the generalized Watson sums, which arise in the theory 
of magnetism and other order-disorder phenomena, in analytical form. These sums 
must be evaluated in order to calculate the thermodynamic quantities of interest as a 
function of temperature. Previous calculations of the generalized Watson sums have 
been mostly based on series expansions and are therefore valid only over limited tem- 
perature regions. The purpose here is to describe a method of evaluating the Watson 
sums without using series expansions or extensive computer summations. The third 
objective is to propose a new type of decoupling. This report shows that it is theoreti- 
cally possible to have different orders of transitions for the anisotropic Heisenberg 
model. 

The simplest way of treating correlations between atoms in a ferromagnetic crystal 
is to consider only one atom and replace its interaction with the rest of the crystal by an 
effective field. This theory, called the molecular field theory, was first advanced by 
Weiss (ref. 1), who assumed the effective field to be proportional to the magnetization. 
This theory gave a satisfactory qualitative picture of the main features of a ferromagnet, 
namely, the existence of a transition temperature, specific heat anomaly, and a sponta- 
neous magnetization. It had weak points, however, such as its failure to predict the 
correct behavior of these quantities at very low temperatures and its underestimation of 
the transition temperature. It also failed to explain the origin of the molecular field. 

Not until the advent of quantum mechanics could the Weiss molecular field be justi- 
fied. Heisenberg (ref. 2) showed that this field is a result of a quantum-mechanical ex- 
change interaction. The exchange interaction is a consequence of the Pauli principle that 
requires the wave function to be antisymmetric with an exchange of space and spin coor- 
dinates. This exchange interaction can be interpreted naively as an interaction which 
tends to orient the spin angular momenta and hence affect the magnetic moments of the 
atoms. 

Dirac (ref. 3) formulated a Hamiltonian which displays the effect of coupling of spin 
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angular momenta of localized particles. This Hamiltonian is often referred to as the 
Heisenberg model, the mathematical statement of which is 

H = -Vj..S. • S . (1) 

H 1 1 

ij 

where J . . is the exchange interaction between the spins on sites i and j and the 

1J fK 

summation is over all sites. The vector spin of the i particle is 

The exchange interaction is positive for ferromagnetism and in practice is regarded 
as a phenomenological constant, which is determined by fitting the theory with a set of 
experimental results. 

In principle, all the thermodynamic and magnetic properties of a crystal that con- 
tains N lattice points can be determined by constructing the partition function for the 
Hamiltonian, and taking the appropriate derivatives. In practice, it has not been pos- 
sible to solve this problem only for certain limiting cases. The only solutions that have 
been obtained for the Heisenberg model are for very high or very low temperatures. 
These solutions have been derived from series expansions of the free energy in powers 
of T (the spin wave solution) valid for T « T and in powers of 1/T for T » T . 
Each of these solutions is only valid far from the transition temperature and therefore 
various approximations have been proposed to study the behavior at intermediate tem- 
peratures. 

The existence of a phase transition in a system with long-range coupling is a ques- 
tion for statistical mechanics. At the present this question cannot be answered mathe- 
matically for the full exchange interaction. But the known (high and low temperature) 
solutions will provide both physical guides of the choice of approximation to be applied 
and the criteria against which to judge these approximations. 

Perhaps the most powerful technique in the study of ferromagnetism is the method 
of temperature-dependent Green’s functions (refs. 4 and 5). With this method the 
thermodynamics as well as the correlations of the spins can be analyzed. Green's func- 
tion theory when applied to ferromagnetism agrees with the spin wave theory at low tem- 
peratures and with the 1/T expansion at high temperatures. 

Some of the complications of the quantum mechanical theory of ferromagnetism are 
avoided if the Ising (ref. 6) model is used rather than the Heisenberg model. If, in the 
Heisenberg model, one assumes that the instantaneous values of the neighboring spins 
may be replaced by their time averages, then one obtains an Ising model. The Ising 
model may, however, be considered the limiting case of the Heisenberg model, in which 
the exchange forces are extremely anisotropic in the z direction, and thus the Hamil- 
tonian is just 
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(2) 


H'-ZVi*? 

ij 


The suppression of the x and y components of the spin in this model produces 
consequences that must be treated with caution. This model is incorrect at low temper- 
atures because it excludes spin waves. The absence of spin waves or low-temperature 
elementary excitations occurs because the quadratic x and y dependence of the energy 
is neglected. The advantages of the Ising model is that calculations can often be carried 
out with rigor and it can be solved exactly in one and two (ref. 7) dimensions. 

The exact solution of the two-dimensional Ising model shows a spontaneous magne- 
tization proving that the system is ferromagnetic. In contrast, the two-dimensional 
Heisenberg model has no ferromagnetism. In three dimensions, by the use of various 
approximations (refs. 8 to 10), the thermodynamics of these models can be estimated. 
The results of these approximations show that both models appear to exhibit a ferro- 
magnetic behavior. 

Besides the contrasting behavior pattern of the Ising and Heisenberg models in two 
dimensions, there- are other apparent differences in three dimensions. Thermodynamic 
quantities like magnetization and susceptibility behave differently below the transition 
temperature. The only place where both models show similar behavior is in the para- 
magnetic state. 

Because of the dissimilar behavior of these models, Dalton and Wood (ref. 11) have 
considered an intermediate model. They have studied the critical behavior of this aniso- 
tropic model in two and three dimensions for intermediate values of an anisotropy param- 
eter ? 7 , where 0 < rj s 1, by the Green function method and series expansion technique. 
The Ising and isotropic Heisenberg models are special cases with values of rj equal to 
0 and 1, respectively. The Dalton and Wood Hamiltonian is 


H = 


ZZ J ij s f s f + 




+ s r s J 


(3) 


where 


S* = S* ± iS7 (4) 

In a physical sense, r/ can be interpreted as a strength parameter for an anisotropic 
exchange interaction. It is usually so interpreted for the case of ferromagnetic crystals 
where crystal field effects give rise to magnetic anisotropy. The magnetic anisotropy 
can have significant effects on the correlation and long-range order properties of the 
system. 
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The exact equations of motion for any order of the Grqen's functions for the aniso- 
tropic Heisenberg model involve higher order Green's functions. These higher order 
Green’s functions thus necessitate a decoupling approximation in order to obtain a solv- 
able closed system of equation. From these equations the correlation function and hence 
the thermodynamics of the system can be exposed. Various types of phase transitions 
are predicted, depending on the form of the decoupling parameter. This will involve a 
solution for the magnetization that will be valid for all temperature regions of interest 
without the use of any separate approximations for high or low temperature regions. 


REVIEW OF MOLECULAR FIELD THEORY OF FERROMAGNETISM 

A qualitative understanding of the Weiss molecular field theory can be obtained by 
examining the Ising Hamiltonian of equation (2). This equation may be considered to be 
the statement that a spin S z feels an effective field 

= < 5 > 
j 

where the sum is over all nearest neighbor sites of the site i, such that the Hamiltonian 
is 


H = 



It is assumed here that the exchange energy between nearest neighbor electrons of 
the same atom is constant and that all electrons have the same exchange energy. Only 
nearest neighbor interactions are included. The molecular field approximation consists 
of replacing equation (5) by 


<^ i >=Jz(S ] Z ) (6) 

where z is the number of nearest neighbors. The process of replacing the spin oper- 
ator S z by its statistical average /s z \ is the essential approximation of this method. 

1 ' ■/ /z\z 

The replacement is equivalent to neglecting fluctuations of J , if is up, its neigh- 
bors will have more than average predilection for being up. The problem is then reduced 
essentially to treating a paramagnetic gas of noninteracting spins. 

As all the magnetic atoms are identical and equivalent, ^S z ) is related to the total 
magnetic moment of the crystal by 
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(7) 


M = N(S Z ) 

where N is the number of atoms and M is the magnetization. Equation (5) can then be 
written in the form 


= — M 
1 N 


(8) 


The thermodynamic quantities can be determined directly from the partition func- 
tion. It is a straightforward matter to calculate the average z component of spin 


S 1 - 


sf exp (.#<*,> 2 s ? 


exp (#<#!> Y, S i 


where 


1 

KT 


Writing out the trace for spin 1/2 and setting ^S z ^ = (s z > yields 


<S Z >=I 


eXP \2 ^ V " 6XP \ 2 ^ 


exp + exp <JSf 4 > 


or 


< S z ) =| tanh^— ^Jz(S z )j 


(9) 


using (,#.) = J(S Z ) . Comparison with various experimental data gives a surprisingly 
good qualitative picture of the thermodynamic properties such as the temperature depen- 
dence of the magnetization, susceptibility, and specific heat. However, the molecular 
field theory does not describe accurately the exact behavior of a spin system near abso- 
lute zero or the Curie point. 

There have been many extensions of and improvements to the Weiss theory for both 
ferromagnets and antiferromagnets. In the molecular-field theory approach, the ex- 
change interactions in the crystal are replaced by an effective field, so that certain 
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interactions arising from S * S are lost. If the exchange interaction were to be consid- 
ered, the problem would become unmanageable. A compromise has been sought where 
a small section of the crystal is treated; the exchange interaction in this region can be 
solved exactly. The effect of the remainder of the crystal was then taken to be an effec- 
tive field. The first detailed treatment of the exchange-coupled spin pair in an effective 
field was given by Oguchi (ref. 12). He considered the effective field to be a molecular 
field which is proportional to the average magnetization. This method was an improve- 
ment over the Weiss field theory in describing short-range order above the critical 
point, and it gave a nonvanishing specific heat. On the other hand, it failed for the tem- 
perature dependence of magnetization near absolute zero and the transition point. The 
Bethe-Peierls-Weiss method showed some improvements over the Oguchi technique as 
far as short-range order, but it could not be used for temperatures much below the 
Curie point because of the expansions used in the effective field. Thus, one of the major 
flaws in these techniques has been the use of expansions that either did not converge or 
gave very inaccurate results. The shortcomings of the effective field theories at low 
temperatures have been corrected by the spin-wave theories. 

The effective field theories are also inaccurate at high temperatures where local 
correlations are important. 


GREEN'S FUNCTION THEORY OF ANISOTROPIC HEISENBERG FERROMAGNET 

In this section, the Green's function theory is applied to the anisotropic Heisenberg 
model, using some of the decoupling approximations previously applied to the isotropic 
Heisenberg model. The algebra of the Green's functions is presented in appendix A. 
The Hamiltonian for the anisotropic model is 


where 


H = J„S?S? + 

i i 


i] i J 


H (sts~ + s:s +N ) 

9 V 1 J iJ/ 


( 10 ) 


S± = S* ± iSj (11) 

77 is the anisotropy parameter, J- j is the exchange interaction between spins on sites i 
and j, and the sum is carried over all sites in the crystal. The exchange interaction is 
assumed to be a function only of the distance between sites. The self- exchange terms 
such as or vanish. 

The anisotropic parameter rj separates the interaction into an isotropic and aniso- 
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tropic part and, in effect, takes into account any crystalline field or lattice distortions. 

Since the thermodynamic quantity of interest is the magnetization of the system, one 
is interested in evaluating correlation functions of the form and hence of the 

retarded Green's function S m})* The reason for this choice be apparent in 

the section on magnetization. 

Using the Hamiltonian in equation (10), together with the known spin commutation 
relations, results in 



One can write the equation of motion (A29) in the form 


E 



a - ( s i) 


5 gm + 2 y J 

w" V ** 


s ? s r s m 


» - "{< 


slsi 
g f 



(12) 


(13) 


For a system of interacting spins, the second and third terms on the right hand side 
of equation (13) contain Green's functions with three operators and are called higher or- 
der Green's functions. The equation of motion for one of these higher order Green's 
functions involves the next higher order Green's function. In this way an infinite hierar- 
chy of equations is generated. Thus the exact treatment of the equations of motion in- 
volves the solution of an infinite set of coupled equations for an infinite number of 
Green's functions. To obtain a closed solution for a Green's function, this hierarchy of 
equations must be truncated at some point. The procedure usually adopted is to make a 
so-called decoupling approximation in which the higher order Green's functions on the 
right hand side of equation (13) are expressed in terms of lower order Green's functions. 
By doing this, equation (13) can be explicitly solved for 


Random -Phase Approximation 

One of the most commonly used decoupling schemes is the random-phase approxima- 
tion (RPA; ref. 13). This technique ignores any fluctuations of and the operator 

S* is replaced by its average value 

o 

(i s l s t ' S m)) IT ( s l) (( S l ; S m» < 14 > 
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A physical justification of this decoupling procedure is plausible at low temperatures 
where the spins are nearly completely alined. It is here that any deviation of the z 
component of the spin at any one lattice site is very small; therefore, neglecting the 
fluctuation is valid for a first approximation. Although this approximation yields quali- 
tatively correct results for the magnetization, it overestimates the transition point. 

Because of translational symmetry of the lattice, the subscript can be dropped in 
(S z ) . On substitution of equation (14) into equation (13), 

[E - 2(SV(0)]A = i <S Z >« - 2 tj(S z ) (15) 

n f 


results where, for brevity, 



(16) 


and 


J(k) =^Jgf exp[ik • (f - g)] 
f 

The translational invariance dictates the use of Fourier transformation with respect 
to the inverse lattice, that is, 


A gm^ ^expfik. (g - m)] A g. 
k 


(18) 


where N is the total number of sites in the lattice and the sums are restricted to the 
first Brillouin zone. It is shown in appendix A that the Green's function for this model 
can be written in the following form 


A k = ~ (S Z )[ E - E(k)]" 1 

77 


where 


E(k) = 2 (S z ) J(0) 


1 


7? J(k) 

J(0L 


(19) 


( 20 ) 


the elementary excitation energy. It should be noted in contrast to the isotropic case 
that the ground-state energy E(0) has a gap of magnitude 2(1-?7)(S Z ) J(0). 


9 



Symmetric Decoupling 


Callen (ref. 14), using a heuristic argument based on physical grounds, suggested 
a decoupling that attempts to account for the spin correlation on one site with the spin 
on another. He proposed a symmetric method of decoupling the Green's function 

(( S ?g S fl S m)) 


such that 


( S ?g) (Ki S m)) + ( S g S g) <Ki S m)) * ( S ??) ({*&)) < 21 » 


s~s + sts ,/ 
g g fj m/I g*f 

where the last term vanishes because it has no diagonal elements in the repre- 

sentation of the total z component of spin. The symmetric decoupling method is simi- 
lar to the Wick-Bloch-de Dominicus theorem for statistical averages (ref. 15). For 
spin 1/2, S z can be written exactly as 


S z = S - S~S + , IS = - 
g gg’l 2) 


( 22 ) 


Furthermore, for any value of spin, from equation (12), 


S Z = - fes: - &~S 

g 2 \ g g g g, 


(23) 


Multiplying equation (22) by a, equation (23) by (1 - a), and adding result in 


“c+ 


SZ= aS + 


(1 - or)SX (1 + a)S~S 


g g 


g g 


(24) 


Substituting equation (24) into equation (14) and using the symmetric decoupling 
scheme (eq. (21)), one finds 


(( 


S Z S/ ; 
g f ’ 





m 


)) 


- (a 


S St 
g f 



g 



(25) 


Here, if 0=1, the decoupling is on the basis of identity (22) and is strictly valid only 
for S = 1/2. If o = 0, the decoupling is on the basis of identity (23), which holds for 
arbitrary spin. Actually equation (22) is valid for arbitrary spin, also, if it operates 
only on spin states that have S = ±S, in effect neglecting deviations from spin-up and 

B 
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spin-down. Thus, the Callen decoupling scheme with a suitable choice of a may give a 
correction to the random-phase approximation by introducing deviations from spin-up 
and spin-down. A physical criterion is required to select the proper value of a. 

The operator S”S* in equation (22) represents the deviation of S z from S. It 
therefore seems plausible to use equation (23) when (S z ) ~ S, that is, at low tempera- 
tures. Thus a = 1 corresponds to a low-temperature decoupling scheme. It should be 
noted that the random-phase approximation neglects these terms. 

On the other hand, as seen from equation (22), the operator 1/2(S + S” - S"S + ) rep- 
resents the deviation of S z from zero, and it is therefore reasonable to use this equa- 
tion when S z = 0. Thus, Callen chose a = 0 to represent a high-temperature decoupling 
scheme. 

For S = 1/2 all the preceding requirements are satisfied by the choice 

a = 2{S Z ) (26) 

When a = 0, one obtains the same equations of motion as for the random-phase approxi- 
mation. However, it should be emphasized again that this is not the only possible choice 
that can satisfy the physical boundary conditions. In other words, a is not unique. In- 
serting equation (25) into the equation of motion (eq. (13)) gives 

EA = £2 + 2(S Z > (A - ti A (m ) + 2a XW gm - VW (2,) 

71 f f 


where again for brevity, 


and 



(28) 


(29) 


When it is assumed that the lattice is translationally invariant, a Fourier transform 
with respect to the reciprocal lattice can be used to obtain the following result (see ap- 
pendix B): 


A - < sZ > 

k tt[e - E(k)] 

where, here, for the Callen decoupling scheme 


(30) 
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and 

J (k) = Yj J gf ex P[ _i ^ • (S- f)] (32) 

f 

The Curie temperature is defined as that temperature at which (S z ) vanishes. 

Then from equation (26) this is the temperature at which a - 0. However, the form for 
a is by no means unique. In this report the decoupling parameter a will take on other 
forms which also meet the physical criteria. It will turn out that, not only is the Curie 
temperature sensitive to this parameter, but so is the type of phase transition. 


MODIFICATIONS OF THE DECOUPLING PARAMETERS 

As mentioned previously, for a second-order transition the long-range order pa- 
rameter (magnetization for the anisotropic ferromagnet) changes continuously from the 
ordered to the disordered state but has a discontinuity in its first temperature deriva- 
tive. In contrast, for a first-order transition, the long-range order parameter itself 
changes discontinuously at the transition point. Other types of transitions appear in 
which the long-range order decreases smoothly with a discontinuity in some higher order 
derivative. There is also the possibility of a transition with no discontinuities in any 
derivative. Such a process will be called a washed-out transition which is characterized 
by a tail in going from the normal to the condensed state. The washed-out behavior is 
believed to occur for the Ising model of solid orthohydrogen (refs. 16 and 17). 

It will be shown that these transitions appear to be quite sensitive to the choice of a 
and /3. In the case of the Callen decoupling, certain criteria were used which imply that 
a must be equal to 1 at T = 0 and must approach zero at T = T . Thus, to meet these 
conditions, the calculations were based on a = 2(S Z ) for spin 1/2. Copeland and 
Gersch (ref. 18) changed the parameter to the form 

a=(2<S z >) P (33) 

where p is any power. Applying this parameter to the isotropic Heisenberg model, they 
reached the conclusion that this form with p = 3 gave better results near the Curie point 
than p = 1 and that the result still agreed with spin-wave theory at low temperatures 
and 1/T expansions at high temperatures. They also showed, with the use of a com- 
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puter, that some powers of p gave a first-order transition and that others were ruled 
out because they do not produce a transition. 

The anisotropic model is applicable not only to ferromagnetism but also to a host of 
other physical problems that involve phase transitions. One would like not only to in- 
vestigate Callen’s parameter, with or without various powers, but also to modify the de- 
coupling scheme such that the critical point values can be improved as well as to yield 
other types of transitions. 

Since there appears a logarithmic variation of the long-range order, a suitable 
modification which is used in this report has the form 

« = (2 (S Z )) P - (2<S Z >) P log(2<S z >) P (34) 

When (S z ) = 1/2, a = 1 at low temperatures, and when (S z ) =0, a = 0 above the 
critical point. This choice also meets Callen’s physical criteria. The results of these 
modifications will be discussed in the next two sections. 

The purpose in modifying the decoupling approximation is to determine what types of 
phase transition can be present. Conventionally, the ferromagnet is believed to undergo 
a second- order phase transition from the ordered ferromagnetic phase to the disordered 
paramagnetic phase at a critical temperature called the Curie temperature. However, 
because of intrinsic difficulties in measuring the thermodynamic anomalies at the criti- 
cal point, the experimentalists are not able to determine what type of phase is present. 
The Heisenberg model has never been exactly solved, and thus there is no rigorous 
proof that the Heisenberg model predicts a phase transition. However, with the approxi- 
mations discussed above the Heisenberg model does exhibit a phase transition. Thus, it 
is of interest to study, within these approximations, at what temperature and in what 
manner the spontaneous magnetization goes to zero. 


SPONTANEOUS MAGNETIZATION AND CRITICAL BEHAVIOR: 
RANDOM-PHASE APPROXIMATION 

Spontaneous magnetization means the presence of a magnetic moment in the absence 
of an applied magnetic field. If there were such a field present, some of the spins in a 
paramagnetic system would line up parallel to it, and there would be long-range order. 
The average z component of the spin can be considered to be a measure of this long- 
range order. The thermal average (S z ) is proportional to the magnetization. 

In order to calculate the spontaneous magnetization (S z ) , one makes use of the 
commutation relations of equation (12) as well as 
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s • s = s‘s + + s z + (s z ) 2 


(35) 


which fixes the magnitude of the spin. Thus the average value of the spin correlation is 

(S"S + > = S(S + 1) - <S Z > - ((s z ) 2 ) (36) 

For spin 1/2 equation (36) reduces to 

(S"S + ) = 1 - <S Z > (37) 

2 

The magnetization for the anisotropic Heisenberg ferromagnet is calculated explicitly 
using the decoupling schemes described previously. 

The calculation of the correlation function follows directly from equation (A30) in the 
limit t = t'. Using the 6 function representation, 


5(x) = J- (_L_ 
27ri \x - ie 

lim 

e—0 


x + ie 


it follows that 


where 


<f> 


(S"S + > = 2(S Z > 


= 1 V 

N / y exp[E(k)/3] - 1 


where E(k) is given by equation (20). 

Using equations (37) and (39) one finds that 


<S'S + > = 2(S Z )<f> = - - <S Z ) 


The expression for (S ) then becomes 


<S Z ) = 2 


1 + 2 $ 


(38) 


(39) 


(40) 


(41) 


(42) 
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And it is easily shown that 


$ = 



cothffi 

2 


1 

2 


(43) 


The magnetization as a function of temperature is evaluated from the sum in equa- 
tion (43) over all values of k in the first Brillouin zone of the appropriate lattice. Ex- 
cept in the low- and high- temperature limits numerical methods are generally used. 

But in general, numerical solutions are difficult. It is a purpose of this report to show 
that analytical solutions of equation (43) are possible. 

Consider only crystals with cubic symmetry. For these lattices replace the sum in 
equation (39) by an integral 


1_ 

N 



v 

(27 r) 3 


/// dir 


where v is the volume per site: 
for the simple-cubic (sc) lattices 


for body-centered-cubic (bcc) lattices 



a 


3 


for face-centered-cubic (fee) lattices 


v = 


1 

4 


a 


3 


(44) 


and the integration is over the first brillouin zone of the appropriate lattice. For the 
cubic lattices, equation (44) can be stated in the form 


1_ 

N 



j- r r r * d y dz 

3 * / 0 0 y 0 


(45) 


where for sc lattices 
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(46) 


x = k x a 


y = V 


z = k z a 


for bcc lattices 


x = - k^a 
2 x 


y = 2 k y a 


z = - k a 
2 z 


and for fee lattices 


x = - k a 
2 ^ 


y 4 v 


COth(77t) = + — 

7 rt 


z = - k a 
2 z 


QO 

mV _ l_ 

Tt 77 / ^ R 2 + t 2 


R=1 


one may write equation (43) in the form 



+ I 


2 


where 


I 


1 


VT 

(2tt) 3 2(S z ) 



1 - ??J(k) 
J(0) 


(47) 


(48) 


(49) 


(50) 


(51) 



Here r(R) = R 7 re/(S Z ); r = kT/J(0), the reduced temperature; and cc indicates the 
complex conjugate. 


Body-Centered-Cubic Lattice 

Consider first the case of a bcc lattice, where (from eq. (32)) J(k) can be written as 
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( 53 ) 


r 


Jfk) k x a k v a k z a 

—Li = cos cos — L_ cos 

J(0) 2 2 2 


From equation (47) it follows that equations (51) and (52) have the form 


Ij(bcc) = 


T 

(2tt) 3 (S Z ) 



dx dy dz 

T] cos x cos y cos z 



dx dy dz 

17 cos x cos y cos z 


(54) 


(55) 


The integrals appearing in equations (54) and (55) are evaluated in appendix C. The 
summation in equation (55) is evaluated by the use of Laplace transforms, and this 
method is shown in appendix D. The results are 


I, (bcc) = — - — 
2<S Z > 


— K(k) 

71 


2 


(56) 


I (bcc) = — L_ 

1 (2Q) 


Q coth(Q) - 



17 2 csch 3 (Q)coth(Q) 


(57) 


where K(k) is a complete elliptic integral of the first kind, 



(58) 


and Q = (S z )/t. On substitution of equations (56) and (57) into equation (50), the mag- 
netization is obtained from equation (42). 

No difficulty arises in the limiting case where rj = 0, nor does it have to be consid- 
ered separately as is the case when low- and high-temperature expansions (ref. 21) are 
used. In fact, by setting rj = 0 in equation (42), a result equivalent to the familiar Ising 
result for magnetization (eq. (9)) is obtained, 


<S Z > 


— tanh Q 
2 


(59) 
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Face-Centered -Cubic Lattice 


The treatment for the fee case is analogous to the bee calculation. The results are 


Ij(fcc) 



(60) 


Ig(fcc) 


Q 2 (S z )g 


gQ coth(gQ) + ^ Q 2 csch 2 (gQ) - — + 
3 3g 



where 


g. 2 Q 3 

18 


esen (gQ)coth(gQ) 


(61) 


,2 _ 4??(3 + rj) cos z (x) 

1 = 9 

(3 + v r 


(62) 


g = 


3 + V 


(63) 


and 


^ (k) _ cos COS — + COS - 2 — COS —5— + COS —5— cos x 

J(0) 


.a k..a 

cos - 
2 2 


k„a ka 

Cl 

2 2 


k a k„a 

Z 
2 


(64) 


Simple-Cubic Lattice 


The results for the simple-cubic lattice are obtained in a similar fashion: 

KU 2 )dx 


^(sc) = 


2 0 

77 Q 


/' 


3-7? cos(x) 


(65) 
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Magnetization, (s z ) 


Anisotropy Anisotropy 

parameter, parameter, 



0 .1 .2 .3 .4 .5 0 .1 .2 .3 


Reduced Curie temperature, t c 

(a) For body-centered -cubic lattice. (b) For face-centered-cubic lattice. 

Anisotropy 

parameter, 



0 .1 .2 .3 .4 .5 


Reduced Curie temperature, x c 
(c) For simple cubic lattice. 

Figure 1. - Magnetization as function of reduced Curie temperature. 


I 2 (sc) = 



coth(pQ) + 


3 


cscti (pQ)coth(pQ) 


dx" 


- J_ (9 - „ 2 )' 1/2 1 jfa - r?)' Vi 

2Q 2 


where 



2 77 

3 - ?? cos(x) 


(67) 


3-77 cos(x) 
3 


( 68 ) 


and 


= [cos(k a) + cos(k a) + cos(k a)l (69) 

J(0) L x y z J 

Figure 1 shows the magnetization (S z ) for each of the cubic lattices plotted as a 
function of the reduced temperature t. Numerical values were obtained from equa- 
tion (42) for several values of 77 . 

Critical Behavior in Random -Phase Approximation 

The random-phase approximation predicts that the spontaneous magnetization shows 
a second-order phase transition independent of the anisotropy. The effect of the aniso- 
tropy parameter is to shift the Curie temperature. The Curie temperature decreases 
smoothly as 77 increases. Physically, this decrease is due to the energy it takes to 
reverse a spin. For the Ising model, since there is only one component present, an 
energy gap of magnitude 2J(0) is needed to create an excitation. In the isotropic model, 
the gap vanished because of spin waves. These physical (intuitive) arguments are sup- 
ported by the numerical calculations presented in this report. 

The Curie temperature is easily obtained from the derived results. For example, 
consider the bcc lattice. By substituting equations (56) and (57) into equation (42), one 
finds 
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1 


2 r = 


K ( k ) 


2 + Q coth(Q) - 1 + 77 ^ csch 2 ( Q)coth(Q) _ 77 ^ 

8 8 


(70) 


where 


and 


T = 


kT 

J(0) 


Q = 


itl 


Z o o 

As S approaches zero, terms such as Q coth(Q) and Q csch (Q)coth(Q) become 
indeterminate. Use of L'Hospital’s rule yields 

lim Q coth(Q) = 1 (71) 

Q-0 


lim Q 2 coth(Q)csch 2 (Q) = 1 (72) 

Q-0 


and 


27 c = 


2K(k) 


-2 


(73) 


Thus, the Curie temperature can be determined without the use of any expansions. The 
same procedure can be applied to the other cubic lattices. 

In table I the values of the reduced Curie temperature r c are given for various 
values of 77 for the three cubic lattices. Also shown are results for r c calculated from 
other approximations. The results for 77 = 1 agree with the Curie temperature of a 
Heisenberg ferromagnet calculated in reference 13 which also uses the random-phase 
approximation. 

The preceding calculations were limited to positive values of 77 . However, for 
these calculations the thermodynamic quantities are independent of the sign of 77 . Per- 
haps, if next nearest neighbors were taken into account, different results would be ob- 
tained. 

Of particular interest is the behavior of the magnetization as 77 varies from the 
Ising to the Heisenberg limit. No singularity or nonconvergence problems arise for 
0 < 77 < 1 in the analytical expressions derived, which are valid for the total tempera- 
ture range. Previously, different approximations had to be made for different tempera- 
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TABLE I. - CURIE TEMPERATURES FOR CUBIC LATTICES, NEAREST 


NEIGHBOR INTERACTIONS FOR SPIN 1/2 


Lattice 

Anisotropy 

Brown and 

Tahir- Kheli 

Callen 

This 


parameter, 

Buttinger 

and ter Haar 

(ref. 14) 

report 


V 

(ref. 29) 

(ref. 13) 





Reduced Curie temperature, t q 


Body centered cubic 

0 

r 



0.5000 


.5 

<0. 2987 

0.3588 

0.4625 

.4830 


1.0 

l 



.3588 

Face centered cubic 

0 

r 



0.5000 


.5 

<0.3500 

0.3718 

0.4666 

.4869 


1.0 

l 



.3718 

Simple cubic 

0 

r 



0.5000 


.5 

< 0. 3166 

0.3296 

0.4500 

.4776 


1.0 

l 



.3297 


ture regions of interest. In fact, as shown in reference 11 77 « 1 produces series 
convergence problems. 

The data in table I shows that the reduced Curie temperatures in the present calcu- 
lation are lower than the results Callen (ref. 14) obtained using a symmetric decoupling 
scheme and are slightly larger than the Brown and Luttinger (ref. 19) results, which are 
based on extrapolation of the high -temperature series expansion. The low -temperature 
expansions for spin 1/2 tend to be in agreement with Brown and Luttinger. The Callen 
case will be discussed in the next section. The anisotropic Heisenberg ferromagnet has 
been described in references 20 and 21. 

At 77 = 0, t c is given by the Ising model values which, as shown in the figures 
1 to 3 and table I, are independent of the lattice structure. 

SPONTANEOUS MAGNETIZATION AND CRITICAL BEHAVIOR IN THE 
SYMMETRIC DECOUPLING APPROXIMATION 

For the symmetric or Callen decoupling approximation one finds that 
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and 


J ( k ) =Y J_f exp[-ik • (g - f)] 


(75) 


By use of equations (A30) and (38), an elementary integration yields the following rela- 
tion for the correlation function: 


* k , =2<S Z >$ 


where 


(76) 


Y i =Y«V 

N Z^i exp[/3E(k)] - 1 L_j K 


(77) 


Thus equation (74) can be replaced by 


E(k) = 2<S Z > j[j(0) - 7?J(k)] + — Yj [7 J (k) - J(k’ - k)]</>J 

Ik* ' 


(78) 


For the nearest neighbor approximation and constant J, 


J(k) = J Y ex P(ik ' 5) (79) 

6 

where 5 is a nearest neighbor vector. It is shown in appendix E that the term 

[j(k’) - J(k’ - k )]<p k , (80) 

which appeared in equation (74), can be replaced by 


Hence, 


[j(0)] _1 |[j(0) - J(k)Y J (k’V^| 

(81) 

E(k) - 2<S Z )J(0) 1 + 2ar]f " J ®. 

(82) 

j (o) (n + 2 off) 
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I 



where 


f 




J(k’) 

J(0) 


<P k * 


(83) 


The definition of E(k) as given by equation (82) together with equation (79) are the 
basic equations of the theory; (p-^, is determined by the requirement that these equations 
be self-consistent. 

The calculation of the spontaneous magnetization (S z ) for the Callen decoupling 
follows in the same manner as in the random -phase approximation. The expression for 
<S Z ) is given by 



1 

2 

1 + 2 $ 


(84) 


where from equations (77), (82), 


and (83), one finds 




1 

exp /3E(k) - 1 



(85) 


where E(k) is given by equation (82). 

Inserting equations (44), (48), and (81) into equations (83) and (85), one obtains 


where 



+ I, 


( 86 ) 
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and 


where 


A = 1 + a(2frj) 

(91) 

B = -q + a (2f) 

(92) 

g = * 

(93) 

A 

f = > 5 + h 

(94) 


From this point on only the bcc lattice will be considered. In appendixes C and D, 
it is shown that the solution of equations (87) to (94) in the nearest neighbor approxima- 
tion is 


Io = 


2<S Z >A L7r 


-K(k) 


(95) 


I/t = 


2(S Z ) A 


Q coth(Q) - 1 + csch2 Q c oth (Q) - G 2 

8 8 


(96) 


I 5 = - 


*6 = 


2QAG 


1 - 


7TK(k) 


Q8A 


Q 2 csch 2 (Q) = 1 


(97) 

(98) 


25 



where K is a complete elliptic integral and 


K 


2 


1 

2 



Vl - G 2 


(99) 


Equations (94) to (99) must now be solved simultaneously for f and E(k). The re- 
sulting value of f is then substituted in equations (84) and (86) to yield the magnetiza- 
tion. 

Figure 2 shows a plot of (S z ) as a function of reduced temperature for several 
values of r) for the bcc lattice using a = 2(S Z ). For 77 = 0, the Callen decoupling ap- 
proximation gives no solutions for a + 0. The only solution in this case is for a = 0. 



Figure 2. - Magnetization as function of reduced tem- 
perature for the body-centered-cubic lattice (Callen 
approximations). Decoupling parameter, a= 2<S Z >. 


TABLE □. - CURIE TEMPERATURES FOR 
BODY-CENTERED-CUBIC LATTICE WITH 
NEAREST NEIGHBOR INTERACTION FOR 
SPIN 1/2 AND a = 2{S Z ) 


Anisotropy 

parameter, 

V 

Callen 

This report 

Curie temperature, r Q , K 

0.01 

Not calculable 

0.5 

.5 

Not calculable 

.4970 

1.0 

0.4625 

.4601 


In table n, the values of kT c /J(0) as calculated from equation (86) are listed for 
several values of V using the Callen parameter a = 2(S Z ) (but not the Callen method, 
which yields only the r/ = 1 result). 


MODIFICATIONS OF CALLEN DECOUPLING PARAMETER 

The formalism for the investigation of the anisotropic Heisenberg model has now 
been developed for an arbitrary a. A choice of the decoupling parameter to be of the 
form proposed by Copeland and Gersch (ref. 18) 
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a 


( 100 ) 


= (<2S Z )) P 

makes it possible to investigate the effect of this parameter on the types of transitions 
and Curie points for various positive powers p. The case p = 1 is just the ordinary 
Callen decoupling which has already been discussed. 

All p ^ 3 give a second-order phase transition at the same Curie point. The Curie 
point value becomes nearly identical to those obtained by the random -phase approxima- 
tion. Figures 3 and 4, which show the reduced -temperature dependence of the magne- 
tization for the bcc lattice, indicate this effect. In contrast to figure 1 for the random - 
phase approximation, the low-temperature behaviors are identical. 



Figure 3. - Magnetization as function of reduced tem- 
perature for a body-centered -cubic lattice. Decoupling 
parameter, a = (2<S Z >). 



Figure 4. - Magnetization as function of reduced tem- 
perature for a body-centered-cubic lattice. Decoupling 
parameter, a = (2<S Z » . 


With o given by equation (100) for 1 < p< 3, one finds a first-order transition for 
the isotropic Heisenberg ferromagnet (?7 = 1) in agreement with the Copeland and Gersch 
computer solution. However, for the anisotropic case, a first-order transition results 
for > 0. 5, as evidenced by double valuedness of (S z ) . In figure 5 the 77 = 1 curve 
shows this first-order transition. This has the consequence that the order of the ferro- 
magnetic transition is not only affected by the decoupling parameters, but also is quite 
sensitive to the value of 77 . 

In table ni the Curie temperature as calculated by equation (84) is listed for several 
values of 77 and p. It should be noted that there is no simple Curie temperature that 
can be obtained from these curves for a first-order phase transition. A transition tem- 
perature for the phase transition must be obtained from a Maxwell construction of the 
free energy for the system. The Curie temperature values in table in for such cases 
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.1 .2 .3 .4 

Reduced Curie temperature, t c 


Figure 5. - Magnetization as function of reduced tem- 
perature for tx)dy-centered -cubic lattice. Decoupling 
parameter, a - (2 S z ) \ 


TABLE III. - REDUCED CURIE TEMPERATURE CALCULATED USING 
a = (2(S 2 ) ) P FOR BODY-CENTERED-CUBIC LATTICE AND 


NEAREST NEIGHBOR INTERACTIONS 


Exponent in 

Anisotropy 

Data from ref. 7 

Data from this study 

Type of phase 
transition 

decoupling 

equation, 

parameter, 

V 

Curie temperature, T^, K 

P 





2 

0.01 

(a) 

0.5000 

2 


.5 

(a) 

.4831 

2 


1.0 

(a) 

.3597 

1 

3 

0.01 

(a) 

0.5000 

2 


.5 

(a) 

.4830 

2 


1.0 

0.3625 

.3589 

2 

10 

0.01 

(a) 

0.5000 

2 


.5 

(a) 

.4830 

2 


1.0 

0.3625 

.3589 

i 2 


a Not calculated. 
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are the intercepts of the (S z ) curve on the temperature axis. 

The formalism has been developed up to this point for an arbitrary functional form 
for a. There is no unique way to determine a within the Green’s function formalism, 
and thus one must rely on experimental observations for this information. 

Structural characteristics of the material have strong effects on the nature of the 
ferromagnetic transition. Inhomogeneity of the concentration, that is, nonuniform dis- 
tribution, of impurities may lead to the appearance of a tail on the curve of the sponta- 
neous magnetization. The ferromagnetic transition can occur, therefore, not at a single 
temperature but over a range of temperatures. Thus, the transition is smeared out. 

One possible reason for this tail is that the inhomogeneity of the concentration leads to 
regions in the material which have different Curie temperatures. In those regions of the 
specimen having a higher Curie temperature, the spontaneous magnetization is still re- 
tained when a part of the ferromagnetic material has gone over into the paramagnetic 
phase (ref. 22). It is suggested, therefore, that the anisotropy parameter can, in effect, 
serve as a measure of inhomogeneity. 

In the section MODIFICATIONS OF THE DECOUPLING PARAMETERS a model was 
developed that had the form 

a = (2<S Z >) P + (2<S Z )) P log(2<S z >) P (101) 

If equation (101) is used in conjunction with equation (74), one finds the following 
results: 

For p = 1 (fig. 6), a washed-out transition is obtained. 

For p = 2 (fig. 7) the isotropic Heisenberg model (?7 = 1.0) gives a first-order 
transition; whereas, for small r/ one obtains a second-order transition. This corre- 
sponds closely to the p = 2 case for the C alien exponential decoupling. 

For p = 3 an important effect comes into play. Previously the exponential Callen 
decoupling parameter (eq. (100)) gave only a second-order phase transition for r\ > 0, 
as well as the same Curie temperature independent of the p. For the present decoupling 
approximation, the order of the phase transition depends on p. For p = 3 (fig. 8) the 
isotropic Heisenberg model (h = 1.0) produces a first-order transition; whereas for 
small ri a second -order transition results. However, for p = 10 (fig. 9) only a 
second-order transition is present. For ij = 0 there is no change in the Ising model 
value for the Curie temperature. Figures 8 and 9 illustrate that various transitions 
occur when the decoupling scheme (eq. (101)) is used. 

Table IV gives the reduced Curie temperatures for this decoupling parameter for the 
bcc lattice. The W means that a washed-out transition occurs. Experimentally, it is 
difficult to distinguish between a higher order and a truly washed-out transition. 



Magnetization, <S Z > 


Anisotropy 

parameter, 



0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

Reduced Curie temperature, r c 


Figure 6. - Magnetization as function of reduced temperature for body-centered-cubic lattice. Decoupling 
parameter, a = 2<S Z > - 2 <S z >ln(2<S z ». 



Figure 7. - Magnetization as function of reduced tern- Figure 8. - Magnetization as function of reduced tem- 
perature for body-centered-cubic lattice. Decoupling peratu re for body-centered-cubic lattice. Decoupling 

parameter, a = (2<S Z » 2 - (2<S Z » 2 ln(2<S z > z . parameter, a = (2<S Z » 3 - (2<S Z » 3 ln(2<S z » 3 . 


r 


■t 
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Figure 9. - Magnetization as function of reduced tem- 
perature for body-centered -cubic lattice. Decoupling 
parameter, a = (2<S Z » 10 - (2<S Z » 10 ln(2<S z » 10 . 

TABLE IV. - REDUCED CURIE TEMPERATURE CALCU- 
LATED a = (2<S Z )) P - (2(S Z )) P log(2(S z )) p FOR 
BODY-CENTERED-CUBIC LATTICE AND THE 


NEAREST NEIGHBOR INTERACTIONS 


Exponent in 
decoupling 
equation, 

P 

Anisotropy 

parameter, 

V 

Reduced Curie 
temperature, 

V 

K 

Order of phase 
transition 

1 

0.01 

0.5000 

2 


.5 


w a 


1.0 


W 

2 

0.01 

0.5000 

2 


.5 

.4841 

2 


1.0 

.3650 

1 

3 

0.01 

0.5000 

2 


.5 

.4830 

2 


1.0 

.3500 

1 

10 

0.01 

0.5000 

2 


.5 

.4830 

2 


1.0 

.3590 

2 


a Washed-out transition. 
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CONCLUSIONS 


From the results obtained by the use of Green's functions and various decoupling 
schemes, certain conclusions can be stated: 

1. A single analytical expression for magnetization, which is valid for all tempera- 
tures, was obtained. The demonstrated technique of evaluating certain Watson sums and 
integrals was employed without requiring various series approximations or extensive 
computer summations. Other thermodynamic quantities, such as specific heat, can be 
calculated from the magnetization. 

2. The random -phase approximation predicts that the spontaneous magnetization 
shows a second -order phase transition at Curie temperature (a discontinuity in the slope 
of the magnetization and in the second derivative of the free energy) independent of the 
anisotropy. The effect of the anisotropy parameter is to shift the critical (Curie) tem- 
perature. The curie temperature decreases monotonically as the anisotropy parameter 
increases. 

3. The Callen approximation for spontaneous magnetization produces self-consistent 
field equations that predict a second -order phase transition. The anisotropy parameter 
shifts the Curie temperature. The values of the reduced Curie temperature are higher 
than those obtained by the random phase approximation. 

4. The exponential decoupling parameter for 1 < p < 3 predicts, by the use of self- 
consistent field equations, a first-order transition for large values of anisotropy param- 
eter h and a second-order phase transition for small values. The type of phase transi- 
tion is sensitive to the value of the anisotropy parameter. For p > 3 (where p is any 
power), a second-order transition is obtained whose reduced Curie temperature is in 
close agreement with the random -phase approximation. The low -temperature dependence 
of magnetization is not the same as that predicted by the random -phase approximation. 

5. The logarithmic decoupling parameter shows what type of transition is sensitive 
to the value of anisotropy parameter. For example, for p = 10 a washed-out transition 
is obtained. But for p = 3 the Heisenberg model gives a first-order transition for large 
values of 77 and a second-order transition for small values. In contrast, p = 2 shows 
the reverse situation. For p = 1 and n = 1, a washed-out transition is obtained. For 

77 < 1, the behavior is that of either a washed-out or higher order phase transition. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 30, 1 970, 

129-02. 
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APPENDIX A 


REVIEW OF GREEN'S FUNCTION 

The Green's function method allows one to focus attention on the quantity of interest 
and to determine from the complete Hamiltonian the exact coupled equations of motion 
for the associated Green's function. In order to solve these equations they must be 
truncated or uncoupled by a suitable decoupling approximation. Although the usual de- 
coupling approximations are not now well understood, this method is a convenient tech- 
nique for the calculation of various thermodynamic quantities as a function of tempera- 
ture. 

In this review, the properties of the temperature -dependent double -time Green's 
function are studied in relation to the work presented in this report. More standard 
treatments of the problems may be found in general review articles and books, such as 
those by Zubarev (ref. 4), Bonch-Bruevich and Tyablikov (ref. 23), Kadanoff and Baym 
(ref. 24), and Perry and Turner (ref. 25). 


Causal Retarded and Advanced Green's Function 

For a system with a time independent Hamiltonian H, the causal retarded and ad- 


vanced Green's function for a pair of operators A and B are defined as follows: 

G c (t, t') = «A(t) ; B(t’))> = -i(TA(t)B(t’)> (Al) 

G R (t, t') = (<A(t) ; B(t')> ) = i0(t - t')([A(t), B(t')]) (A2) 

G A (t, t’) = <<A(t); B(t’))> = i0(t' - t)([A(t), B(t')]) (A3) 

where A(t) is the Heisenberg representation for the operator A, 

A(t) = exp(AHt)A exp(-AHt) (A 4) 

and T is the time ordering operator of Dyson, which is defined in the usual way so that 

TA(t)B(t’) = 6>(t)A(t)B(t') + 0(-t)B(t')A(t) (A5) 

with d(t) as the step function, 
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1 , 


t > 0 


e(t) = 


t < o. 


+i 

[A, B] = AB - crBA: a = 

-1 


(A 6) 


(A7) 


Note that A and B are not restricted to be either Bose or Fermi operators and that the 
choice of a is arbitrary. For the case of Bose operators that are used herein, a = 1 
and the singular angular bracket denotes the average with respect to the canonical density 
matrix of the system 


(O) = Tr^O 


(A 8) 


where 


^ _ exp(-/3H) 

Tr exp(-/3H) 


(A 9) 


and 


f3 


kT 


(A 10) 


Using the equation of motion 


— = [A, H] = AH - HA 


(All) 


and the property of the step function 


dp (t) _ d 
dt dt 



(A 12) 


results in 


^G R (t, t’) - 5(t-t')([ A(t) , B(t’)j) + (/[A(t), H(t)J; B(t’))) (A13) 

and similarly for the causal and advanced Green's function. Note that the Green’s func- 
tion on the right side of the equation can be of higher order than on the left. Writing down 
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the equations of motion for the higher order functions, a set of coupled equations for the 
complete hierarchy of Green's function can be obtained. It is easily shown that the 
Green's function G(t, t’) can be written in the form G(t - t’), and thus depends only on 
the time difference. 

The time correlation function plays an important role in statistical mechanics and is 
usually the quantity of direct physical interest. It is written as 


(A 14) 

which can be calculated from equations (15), (16), and the cyclic property of the trace. 


F AB (t,t') = (A(t)B(t’)) 
F BA (t, t’) = (B(t’)A(t)) 


SPECTRAL REPRESENTATION 


To obtain the spectral representation of the time correlation functions the eigen- 
function m and eigenvalues E of the Hamiltonian must be considered. Thus 

H |m)= E m |m) (A15) 

From equation (A 14), the correlation function can be written 

f BA (ft) = (B(t’)A(t)) 

_ Tr[exp(-/3H)exp(-iHt’)B exp(iHt)A] 

Tr exp(-/3H) 

= Z" 1 Yj ( m l B l n >< n l A |m)exp(-E m /3)exp - [i(E m - E n )(t’ - t)] (A17) 

m , n 


with 

Z =£exp(-/3E m ) (A 17) 

m 

the partition function for the canonical ensemble. Equation (23) can be rewritten in the 
following manner: 
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where 


(B(t')A(t)) = f J(cr)exp[-io;(t - t’)]dcu 

•/-co 


JM = Z" 1 Yj (m|A|n)(n|B|m>exp(-^E m )6(E n - E m - w) (A19) 


In a similar manner 


/* co 

(A(t)B(t’)) = / J(w)exp(o)B)exp[-ico(t - t’)]du> 

•/-CO 


(A 20) 


Equations (A18) and (A20) are the spectral representation for the time correlation func - 
tion, and J(u>) is the spectral density. 

The rest of the discussion will be limited to the retarded Green's function, thus the 
subscript R will be omitted. The Fourier transform of G(t - t’) is 


G(E) = — f G(t - t')exp(iEt)dt 

2tt J-<*> 


(A21) 


Using equations (A2) and (A20) and integrating equation (A21) yields 


G(E) = 


/ICO 

- / [exp(o>/3) - a] ■ 

2ir J (E - a) + ie) 

t/_O0 


(A22) 


where E — 0 . 

Using the symbolic identity 


(w - u>' ± ie) I (w - w') ± in5{u) - w') 


(A23) 


where E — 0 and P denotes Cauch's principal value of the integral, equation (A20) can 
be written alternatively as 


[ exp(fiLo) - a]J(o>) - - [exp(/3E) - cr]J(E) (A24) 

E - co 2 
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It follows from equation (A23) that 


Re G(E) = — 

7 T 



ImGMdcu 

o> - E 


(A25) 


If the function G(E) is known so the spectral intensity J(u>) of equation (A20) can be 
found by taking the difference of the two expressions of equation (A22), then 


G(o)+ie) - G(cu-ie) = 


if 


[exp(/3E) - cr]j(E) 


(cu-E + ie) (w-E-ie)J 


dE 


If the 6 function representation 


(A 2 6) 




6(x) = 

(27ri) 

lim 

e-0 


1 1 

(x - ie) (x + ie) 


is used, then 


(A27) 


G(cj + ie) - G(w - ie) = -i[exp(o>/3) - ct]J(oj) (A28) 

The equation of motion for the Green's function can be obtained by the Fourier transform 
of equation (20) as 


E ((A ; B» =^-(l A > B J) + ((LA, H]; b)) 


(A29) 


Except for the simplest cases, an infinite number of coupled equations are necessary to 
solve for the correlation function (BA) . Thus, these equations must be reduced to a 
finite number by applying an appropriate decoupling scheme. 

The time correlation function is found from 


(B(t’)A(t)) = lim i 
e-0 



G(cu + e) - G(cu - e) exp 


[icu(t - t T )]do) 
exp(/3co) - o 


(A 30 ) 
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APPENDIX B 


FOURIER TRANSFORM OF THE GREEN'S FUNCTION 


Consider the equation of motion for the Green’s function A gm (eq. (27)), which has 
the form 


EA gm = <S*> JB2 ♦ 2<S Z ) £j gf (A gm - VA tm )+2« 2V^gf A g m - 


Defining the Fourier transforms of the Green's function by 




and of the correlation function by 


% = iZ exp E iir ‘ ( ^- 


gf N 4r 

k 


one obtains from equation (Al) 


^exp[ik- (g - m)]A_ = (S z ) 


jexp[ik- (g - mj]A k - ?? exp[ik- (7 - m)]Aj 


+ “ X 2 Yj | 77 ex PQ^* (£"" ™)] ex P "k’ • (g - mjjA^ 

N 2 f k k' ^ 


exp^-ik’ • (7 - gT|exp[ik • (7 - in 3V. A r 

k k 


II i III III l III I ll ll III Mill I l in inia i Hfl 



If equation (B4) is multiplied by 

exp[-ik" • (g - m)J 
g 

and if the following relations are used 

s-Z Z exp C l(ir - £i) - (? - 5 0v E 6 kiT" A k = A k- 

kg k 

Yj 5 gm ex P[- ik " ■ (g - m)] = 1 
g 

^SSE J gf eXp ^‘ (g- m)]exp[-ik". (g- = A, J(0) 

f k g f 


N X22 J gf exp [ ik( * ' m)]exp[-ik”(g - mj\Ag = ^J gf exp[ik(g - f )]a~, 
f k g f 


^ZZZZ^' (? - m)jexp[ik' . (g - f )Jexp^-ik” 
N f k k’ 


• (g - m)]A k ^ k , 


1_ 

N 


ZEZ^E^'-s-^Wkk 

f k k’ 


N 


Z J(ir,) v 

i.' 


k” 


(B5) 

(B6) 

(B7) 

(B8) 

= J(k”)A ? „ 
(B9) 


(BIO) 
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_ N» 

X/X/ X/X/ J gf exp E^‘ m)]exp^ik’ • (f - g)] exp^ik" • (g - m)]A-.^_ 

N f k k’ g 


= ~Y J(k" - k')^ r A p 

xr k k 

k’ 


then 


EA r =iS!l + 

k 


2(S z )A^[j(0) - n J(k)]+ — A~y ??[j(k) - J(k - 

k' 


is obtained by substituting into equation (Bl). Finally, collecting terms of 
equation (30). 


(Bll) 


(B12) 


gives 

k 
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APPENDIX C 


EVALUATION OF INTEGRALS I 
Body-Centered -Cubic Lattice 


1 


It follows from equation (54) that the integral I^(bcc) is written as 

Ij(bcc) = — fff dx dy dz[(l - 77)cos(x)cos(y)cos(z)]^ //2 (Cl) 


Performing the integration on x results in 


Ij(bcc) = — f f dy dx 
1 2 SO SO 

77 


(1 - Tj 2 )cos 2 (y)cos 2 (z) 


1/2 


(C2) 


2 2 

The integrand is expanded in ascending powers of cos (y)cos (z) and then integrated 
term by term to obtain 


oo 

2 /*7 r y*7T 2 2n 2n 

I n (bcc) = — > / / cos(x)cos(y)dy dz 

1 2 / jS§ Jo 

77 

n= 


(C3) 

n= 


where 


= 1 • 3 • 5 . . . (2n - 1) 
2 • 4 • 6 . . . 2n 


(C4) 


where a^ = 1. 

Watson (ref. 26) has shown that the square of the complete elliptic integral of the 
first kind can be written as 


K 2 (k) = 1 7r 2 ^a 2 (2kk’) 2n 


m= 


(C5) 
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where 


/ 

k’ = (l - kj 

Equating equations (C5) and (C3), thus, 

V 2n = (2kk') 2n 

it is easily shown that 


k 2 =I 


/ 

o\l/2l 

H 

H-i 

1 

1 


This gives 


Ij(bcc) = 


2K(k) 

7T 


in agreement with Equation (56). 


Face-Centered -Cubic Lattice 

The fee case can be investigated by the triple integral 

3 nn rn /»7r 

I^(fcc) = /-j III 1 - — (cos(x)cos(y) + cos(y)cos(z) + cos(z)cos(x)j 

11 Jo Jo Jo 


On integrating with respect to x the result is 


IjCfcc) = 


f If (S 

r -n Jo Jo 


- 1/2 


2 \ ' 

a Q + sl 1 cos y + a 5 sin y) dy 


dz 


where 


(C6) 


(C7) 


(C8) 


(C9) 


* dx dy dz 
(CIO) 


(Cll) 
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Let 


~\ 


a »' - 1 


aj = - 2 il±_3) cos(z ) 


a 5 = sin 2 (z) 




Then, 


t = tanf^ 
\2 


2 dt = sec 2 /-jdy = (1 + t 2 )dy 


1 - tan 


m 


cos(y) = 


2/ 1 - t £ 


1 + tan 2 /|) 1 + t £ 


sin(y) = 2 tan^^ 1 + tan 2 /^ = — — 


2 / 1 + t 2 


If these relations are substituted into equation (Cll), then 


Ijffcc) = 2[(„ 2 ,)] ( b 0 )~ 1/2 (a + Bt' 


2 4 \"l/2 

A + r dt dz 


where 


b o = (a o - a l> 


a o + a l 


and 


(C12) 

(C13) 

(C14) 

(C15) 

(C16) 

(C17) 

(C18) 

(C19) 
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(C20) 


B = 


2(a.Q + 2a.g) 


b 


o 


If 


5 


2 





2 



(C21) 


(C22) 


then 


i«cc = [2(^)1 /”/ 


- 1/2 


dt dz 


(C23) 


Once the radical has been expressed in this factored form as the product of sums or 
differences of squares, the integral can be reduced to Jacobian or elliptical form. The 
integral of (C23) is finite everywhere and leads to the integral of the first kind. Thus, 
from elementary tables 


i x ( fcc) = -f- £ ( b o)' 1/2 (^)' 1 K(k i )dz (c24) 

77 7 ? 


where 



4?7 


3 + r) cos^(z) 


(3 + -nY 


Some algebraic manipulation gives 


IjCfcc) 


e_ r K < k i> 

•Vo 3 + ” 


dz 


(C25) 


(C26) 


which is the same as equation (60). 
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Simple-Cubic Lattice 


The integral for the sc lattice in three dimensions can be written as 


I x (sc) = — 


— I II dx dy dz\-— -JJ [cos(x) + cos(y) + cos(z)l 

*Vo Jo Jo 13 


-1 


The same procedure that was outlined in the previous section gives for the fee 


Ij(sc) = — y ^[3 - 77 cos(z)]K(k 2 )dz 


where 


k 2 - 
K 2 " 


2r? 


3-h cos(z)J 


which is of the same form as in equation (65). 


(C27) 


case: 

(C28) 


(C29) 
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APPENDIX D 


SUMMATION AND THE LAPLACE TRANSFORM 


The calculation of the integral I 2 involves the summation of an infinite series. A 
direct approach to the problem of the summation of infinite series has been given by 
several authors (refs. 27 and 28). Some of these techniques will be used and extended to 
evaluate I 2 . 

As an example of the technique used, consider the fee lattice. (The calculation for 
the other two cubic lattices follow in the same manner): the fee lattice integral is writ- 
ten as 



3J(x, y, z) _ CO s(x)cos(y) + cos(y)cos(z) + cos(z)cos(x) (D2) 

J(0) 


The triple integration and the summation exist and are real. However, they can be 
written as the difference of two complex integrals as was done in equation (52). Inte- 
grating on x and y (appendix D) gives 


I 2 (fcc) 



dz 

3 + V + 3iT(R) 


+ cc 


(D3) 


where 


i2 4??[3 + rj cos 2 (z) + 3iT(R)] 

kj - = - - 

[3 + rj + 3iT(R)] 2 


Q = 


r 


(D4) 


(D5) 
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r(R) = — 
Q 


(D6) 


and cc indicates the complex conjugate. 

The complete elliptic integral of the first kind can be represented in series form as 


where 


K(k) = H 

2 



k | < 1 


(D7) 


a = 1 • 3 ■ 5 . . . (2n - 1) 
n 2*4-6. . . 2n 


(D8) 


Only the first two terms of the series in equation (D7) are necessary for four place ac- 
curacy. There is no difficulty in extending it to as many terms as needed. 

Taking the first two terms in equation (D3) and integrating the complex function with 
respect to z, the terms can be rearranged to show that 


where 


I 9 (fcc) = l| + Ig + Ig 


I2 = <1tt) 


00 

I 


R=1 


y 2 + R 2 



2R 



y 2 * * 2 (y 2 + R 2 ) 2 _ 


4R^ 


G = — V = 2 + ^ 

Q G 


(D9) 


(DIO) 


(Dll) 


(D12) 


(D13) 
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The following Laplace transforms are used 


- = if- sin(yt) 

y 


2 w 2 
y + R 


R 


(? 2 * * 2 ) 


if|— [sin(yt) + yt cos(yt)]! 


2 1 2y 


(D14) 


(D15) 


R 


= st< 


(y 2 + R 2 )^ L 8 y 
Equation (D14) can be written as 


3 L 


(1 + y 2 t 2 )sin(yt) - yt cos(yt) 


- = - J sin(yt)[exp(-tR)] 


dt 


2 d 2 y * / 0 

y + R 3 


(D16) 


(D17) 


and similarly for the other terms. 
The summation then becomes 



R=1 



sin(yt) — 

exp(t) - 1 


— — coth(y7r) - — 
2 y 2y 


Using the general form of integration such that 

x2m 


/ 


sin(yt) 


dt 

exp(t) - 1 


(-D 


m 3 


2m 


3y 


2m 


— coth(y7r) - — 
2 2y 



cos(yt) 


dt 

exp(t) - 1 


(-l) m 


Jm+1 

V 

3y 2m+1 


— coth(y7r) - 1 
2 


2y 


(D18) 


(D19) 


(D20) 


results in a closed form solution as was done in SPONTANEOUS MAGNETIZATION AND 
CRITICAL BEHAVIOR: RANDOM-PHASE APPROXIMATION. 
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APPENDIX E 


NEAREST NEIGHBOR INTERACTIONS 


Consider the sum 


K’ = ^ [j(k’) - J(k* - k )]cp~ (El) 

k’ 


which appears in equation (78). For lattice with cubic symmetry and with nearest neigh- 
bor interactions, 


ij 


= J 


if i and j are nearest neighbors, and 


J ij 


= 0 


otherwise. Thus 


J(k) = J ^Texp(ik- 5 ) (E2) 

5 


can be written where 5 is a nearest neighbor vector. Then equation (El) can be written 
as 



exp(-ik- 5)] ^ j exp(ik' 


5)cp r 

k 


(E3) 


By symmetry, the sum over k' is clearly independent of 6, so that it can be replaced 

by 


(lz) I exp(ik • 6 ' ) 
5 ' 


(E4) 


where z is the number of nearest neighbors. Hence, equation (E3) yields 
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K 


' = — ^ [l - exp(-ik 


ci exp(ik' • 6')<p-- t 


6' k’ 


= J 


5 


[l - exp(-ik- 6)] 



J(k') 

Jz cp_r 
k’ 


= — [j(0) - J(k)] 
J(0) L 





since for cubic lattices 


J(k) = j(-k) 


and 


J(0) = Jz 

Equation (82) is obtained by substituting equation (E5) into equation (81). 


(E5) 


(E6) 


(E7) 
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